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Abstract 

In this paper, spin-dependent transport through a spin diode com- 
posed of a quantum dot coupled to a normal metal and a ferromagnetic 
lead is studied. The current polarization and the spin accumulation are 
analyzed using the equations of motion method within the nonequilibrium 
Green's function formalism. We present a suitable method for computing 
the Green's function without carrying out any self-consistent calculation. 
Influence of coupling strength and magnetic field on the spin current is 
studied and observed that this device can not work as a spin diode under 
certain conditions. 

1 Introduction 

Spin-dependent transport through a quantum dot (QD) has attracted increasing 
attention in recent years due to development in constructing spin-based devices 
[U [2j [3l m [5]. Study of transport through QDs obtains interesting informa- 
tion about novel physics phenomena such as Kondo effect O El |8], spin and 
Coulomb blockade effects [9l [TOl HH [l2], spin valve effect [13l [M], tunneling 
magnetoresistance [El [16] , zero-bias anomaly [17] and etc. Couphng the QD 
to the different contacts can lead to form different devices for example the spin 
filters constructed by coupling the QD to the normal metals [181 [H] and spin 
diodes constructed by attaching the QD to a normal metal and a ferromagnetic 
lead [ID]- 

Recently, spin diodes have been studied both experimentally and theoreti- 
cally. C. A. Merchant and N. Markovic [20j observed diode-like behavior in a 
carbon nanotube coupled to a ferromagnetic and a nonmagnetic lead. The spin 
diode-like behavior was also reported by Ivan and co-workers [2V in resonant 
tunneUng sandwiched by tunnel barriers with different spin-dependent trans- 
parencies. In addition, spin diode devices were theoretically analyzed in a few 
articles. Souza and co-workers j22, studied a semiconductor QD coupled to a 
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normal metal and a ferromagnetic metal and I. Weymann and co-workers [23j 
studied a spin diode composed of a carbon nanotube coupled to a normal metal 
and a ferromagnetic metal. Both of them used rate equations [231 [13 which 
are valid in the limit Tq << kT that Tq and T are coupling strength and 
temperature, respectively. Here, we use the nonequilibrium Green's function 
formalism (NEGF) [35] also valid in the limit Tq ^ kT to analyze spin diode 
behavior composed of a quantum dot. The Green's function is obtained by 
means of the equations of motion method within the nonequilibrium Green's 
function formalism up to the second order of Hartree-Fock approximation |27| . 
We extract an analytical relation for the electron density without performing 
any self-consistent calculation. The obtained result can be useful for studying 
charge transport through mesoscopic systems by means of equations of motion 
method within Green's function formalism. The influence of spin splitting on 
the performance of the device is investigated. This splitting can be created by 
an external magnetic field or coupling of the dot to a magnetic substrate. In 
addition, the influence of coupling strength on the diode behavior is analyzed. 
Such an examination is impossible by means of the rate equations because it is 
just valid in the limit of very weak coupling. 

The article is organized as follows: the model Hamiltonian and main formulas 
are presented in section 2, in section 3 we present numerical results and discuss 
the difference between results obtained from NEGF and rate equations. In the 
end, some sentences are given as a conclusion. 

2 Description of the model 

Hamiltonian describing a quantum dot coupled to a ferromagnetic and a non- 
magnetic lead is written as 

a,k.a (J 

+ XI \^o.,K'yc\k,o<^a + Tl^k,„d\co,^kA (1) 

where (^^^. ^ica.k.a) creates (annihilates) an electron with momentum fc, spin 
CT in the lead a. = eg + o''^ is the energy level of the dot, A = Iq^ibB 
denotes Zeeman splitting, B is magnetic field and a is equal to 1(— 1) for t (i), 
respectively. d\(df,^ is the creation (annihilation) operator in the dot and 71^ = 
dJ.do- is the occupation operator. IJ is on-site Coulomb interaction strength, 
Ta,k,a describes tunneling between the dot and the lead a and it is also assumed 
that the electron spin is conserved during the tunneling. 

In order to analyze the system, the nonequilibrium Green's function method 
has been used. The retarded Green's function of the dot is GJJ. — —i'd{t — t') < 
{dcr{t), dl.{t')} > that in the steady state it only depends on the time difference 
r = t—t'. Hence, it is better to use its Fourier transform G^(e) =<< o?cr, d^. >>£. 
Using the equations of motion technique for the nonequilibrium Green's func- 
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tion up to the second Hartree-Fock approximation, the spin-dependent Green's 
function of the dot is given by 



6-e. + i(r^(6) + r«(e)) 

<n^ > 



(2) 



where r"(e) = 27r^^ \Ta,k.a\^S{e — ea,k,a) is the couphng strength giving rise 
to the broadening of the dot levels due to tunneling through the left and right 
leads and a stands for the opposite spin a. In the following we use the wide 
band limit i.e. the energy independent broadening r"(e) = F". The electron 
density is given by [17] 

that /"(e) denotes the Fermi distribution function of the lead a with the 
chemical potential fia- Although it seems Eqs.(2,3) should be solved in a self- 
consistent manner to obtain the exact Green's function, we show the Green's 
function and the electron density can be analytically computed. From Eq.([2]) 
the retarded Green's function can be written as follows EHl 



Glie) ^ A^ie) + B^ie) < > (4) 



where 



' - + 4(E„ r;))(6 - - [/ + i(E<. rj)) 

If Eq.Q is substituted into Eq.(I3|), it can be easily shown that the spin- 
dependent density is given by 

< >= n — (6) 

where 

^ f de r^/^(e)-fr^/«(e) 

= — Y^ + TR I^iBAe)) (7b) 

Knowing the electron density Eq.®, the Green's function is obtained from 
Eq.(|3]). Now, we are able to compute the current given by (e = ?i = 1) 

(8) 
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For simulation purpose, we set U — AmeV, kT — 212meV which is large enough 
to guarantee no Kondo effect and = ImeV + a A that A is equal to 1.16 x 
10~'^meV or 0.17AmeV for B — lO^^T and 1.5T, respectively. The coupling 
strength for the metal lead is set r-j- = Fj^ = Fq and for ferromagnetic lead 
^tii) ~ ^o{^ ± p) that p stands for the spin polarization degree of the lead. 
The chemical potential of the left lead (normal metal) is equal to zero and for 
right one, we set — —V, so when the bias is positive the left lead acts as an 
emitter and in the negative bias it acts as a collector. 

3 Simulation results 
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Figure 1: The current polarization as a function of bias in weak magnetic field 
(solid) and strong magnetic field (dashed). The parameters are p = 0.5 and 
Fq = SO^eV. The left and right insets show the QD levels and the chemical 
potentials of the leads in the regions 1 and 2, respectively. Energy difference 
between t and | levels depends on the magnitude of the magnetic field. 

Fig. 1 shows the current polarization ^ = j^"^^ as a function of the bias 
voltage in different magnetic fields. It is observed that in positive bias when 
/x^ < Ea and Ea < fi^ < Ea + U (the dot is singly occupied) the current 
polarization becomes zero because in positive bias the normal metal lead acts 
as an emitter and the spin-up current is equal to spin-down current due to 
F^ = Ff . Note that when the QD is singly occupied, the current polarization 
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depends on (T" — T°) that a denotes the emitter lead [12] • In negative bias, 
the ferromagnetic lead acts as an emitter and the current polarization becomes 
maximum when Sc < < + U because of > so that the spin-up 
current is bigger than the spin-down current. This behavior suggests that this 
system operates as a spin diode in a definite voltage range. We also observe the 
current polarization behaves completely different in the response to weak and 
strong magnetic fields in the regions 1 and 2 shown in the fig. 1. Although our 
results and the results given in j22! are nearly the same in weak B and Fq, they 
are completely different in strong field because of removing degeneracy [30] . 
The left inset can help us to understand the current polarization behavior in 
the region 1 where the ferromagnetic lead acts as the emitter and the dot levels 
are outside the bias window. In the weak magnetic field [B ~ 10^'^T) and 
£^ are nearly degenerate so that the spin-up can be occupied faster than the 
other spin due to > and as a result /-[■ > /j,, therefore, ^ is positive. In 
the strong magnetic field {B = l.ST), the levels of the spin-up and spin-down 
are split so that the is more energetically accessible hence > and as 
a result the current polarization becomes negative. The right inset describes 
the position of the chemical potentials of the leads and the energy levels of the 
dot in the region 2 in which the normal metal lead acts as the emitter. In 
weak magnetic field both levels have the same energy the t electron can leave 
the dot faster due to r:j^ > therefore is bigger than /j, and the current 
polarization is positive. But in the strong magnetic field the \, electron state is 
occupied faster because it is more accessible and as a consequence, /j, > /-f- so 
that the current polarization becomes negative. When \V\ > Sa + U that the 
dot is doubly occupied the current polarization will be constant due to interplay 
between the spin accumulation and the electron-electron interaction. 

Fig. 2 plots the spin current Ispin ^ h ^ versus the bias voltage for 
different coupling strengths. As we expect with increase of Fq the current is 
enhanced because of faster tunneling. We observe when \ea\ < \eV\ < \ea + U\, 
the spin current as well as the current polarization becomes zero just in the 
limit of low or intermediate Fq and with increasing the coupling strength the 
current polarization deviates from zero and as a result the device can not work 
as a spin diode. The dependence of the spin current on Fq is shown in the inset 
in the conditions that the dot is singly occupied. It is observed that NEGF 
and rate equation obtain the same result for Fq < 0.05mey, but the latter one 
could not describe the system well for Fq > O.ObmeV. Indeed, the rate equation 
predicts ^ = for any value of Fq. The unexpected positive and negative spin 
current has been also observed in the regions 1 and 2, respectively. Note that 
the current is measured from the left lead i.e. the current is positive in the 
positive bias. The positive spin current in the region 1 means that /j, is bigger 
than /-[■ because of ej, < £-[•. The negative spin current in the region 2 shows 
that /|. > /-|- due to the same reason. 

Fig. 3 shows the spin accumulation m = n-j- — nj, as a function of the bias 
voltage. It is observed that in the negative bias the spin accumulation is positive 
because in this situation the right lead acts as the emitter and the t electron is 
injected faster and hence the spin-up population is more. In negative bias that 
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Figure 2: The spin current versus bias in different Fq. We set B = 1.5T. Inset 
shows variation of the spin current as function of broadening in V — 9mV. 



the right lead acts as the cohector, the spin accumulation is negative because 
a J, electron has to stay inside the dot for a longer time. Such a behavior was 
recently reported using the rate equations method. Here, there is a significant 
difference between the results presented in the Ref [35] and our results obtained 
using the Green's function method. In p — 1 that the right lead is a half metal 
and in the positive bias when the energy level of the dot is inside the bias 
window, the spin accumulation is equal to -1 (the gray line in the fig. 3b) in the 
weak broadening (Fq = 10 fieV). Indeed, the J, electron injected from the left 
lead can not leave the dot and since the second level + U is outside the bias 
window, the dot contains only a J, electron. But in the other case Fq = 200fieV, 
we observe the spin accumulation is about -0.8. This difference originates due 
to different density of states (DOS) in these couplings. In the weak coupling 
(fig. 3c) the DOS of the t electron has a peak outside the bias window but the 
peak of the DOS of the J, electron is inside the bias window and as a result the 
4, electron only exists inside the dot. In the strong coupling (fig. 3d) the DOS 
is broadened and the DOS of the t electron has a peak inside the bias window 
therefore the fraction of the t electron can exist inside the dot leading to the 
decrease of the spin accumulation. 
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Figure 3: The spin accumulation as a function of polarization in (a) strong 
magnetic field and (b) weak magnetic field. We set p = 0.2(solid), p = 0.4 
(dashed), p = 0.7 (dotted) and p = 1 (dash-dotted). Fq = 200^eV in all plots 
except the gray line in (b) equal to lOfieV. (c) and (d) show DOS (spin-down 
(dashed) and spin- up (solid)) in weak and strong couplings, respectively. 



4 Conclusion 

In this paper, spin-dependent transport through a quantum dot attached to a 
normal metal and a ferromagnetic lead is studied by using the nonequilibrium 
Green's function formalism. We examine under what conditions this device can 
work as a spin diode. We also investigate the infiuence of magnetic field on the 
current polarization. It is observed that coupling strength has an important 
role in the performance of the device. More specifically, in the strong coupling 
this device can not work properly. It is also observed when the dot is coupled 
to a half metal lead, the spin accumulation is equal to -1 in the weak coupling 
but it deviates from -1 in the strong coupling. 
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